EPJ manuscript No. 

(will be inserted by the editor) 



Infrared regularization for spin-1 fields 



o 
o 

(N 

> 
O 



Peter C. Bruns 1 a and Ulf-G. MeiBner 12 b 

1 Helmholtz-Institut fiir Strahlen- und Kernphysik (Theorie), Universitat Bonn, Nufiallee 14-16, D-53115 Bonn, Germany 

2 Forschungszentrum Jiilich, Institut fiir Kernphysik (Theorie), D-52425 Jiilich, Germany 

Received: date / Revised version: date 

Abstract. We extend the method of infrared regularization to spin-1 fields. As applications, we discuss the 
chiral extrapolation of the rho meson mass from lattice QCD data and the pion-rho sigma term. 
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1 Introduction 

Chiral perturbation theory (CHPT) is the effective field 
theory of the Standard Model at low energy [1,2]. It is 
based on the spontaneously broken approximate chiral 
symmetry of QCD. The pions, kaons and the eta can be 
identified with the Goldstone bosons of chiral symmetry 
breaking. Their interactions are weak and vanish in the 
chiral limit of zero quark masses when the energy goes to 
zero. This is a consequence of Goldstone 's theorem and 
allows for a consistent power counting. Consequently, any 
amplitude can be written as sum of terms with increasing 
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powers in external momenta and quark masses, symboli- 
cally 

A = P D f(p/fi,g) , (l.l) 

where / is a function of order one, p collects the small 
parameters, D is the chiral dimension, fi a regularization 
scale (related to the UV divergences in loop graphs) and g 
a collection of coupling constants, the so-called low-energy 
constants (LECs). Weinberg [3] first established this power 
counting, in particular, the expansion in p (the chiral ex- 
pansion) can be mapped onto an expansion in terms of 
tree and loop graphs, with n loop graphs being parametri- 
cally suppressed by powers of p 2n compared to the leading 
trees. The explicit expression for D reads: 



D(A) = V n( n ~ 2 ) + 2L + 2 



(1.2) 



N 



with L the number of Goldstone boson loops and V n a ver- 
tex of order 0(p n ). In essence, this power counting works 
because the pion mass vanishes in the chiral limit and 
thus the only dimensionfull scale in this limit is the pion 
decay constant F v (more precisely, its value in the chi- 
ral limit). Utilizing a symmetry preserving regularization 
scheme like e.g. dimensional regularization leads to homo- 
geneous functions in the small parameters (for more de- 
tails, see the reviews [4-8]). The precise relation between 
this effective field theory (EFT) and the chiral Ward iden- 
tities of QCD was firmly established in Refs. [9, 10]. 

The active degrees of freedom in CHPT are the Gold- 
stone bosons, chirally coupled to external sources. It is, 
however, known since long that vector and axial-vector 
mesons also play an important role in low energy hadron 
physics, as one example we mention the fairly success- 
ful description of the pion charge form factor in terms of 
the (extended) vector dominance approach (for more de- 
tails, see e.g. the reviews [11, 12]). These heavy degrees 
of freedom decouple in the chiral limit and at low en- 
ergy from the Goldstone bosons. Nevertheless, they leave 
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their imprint in the low-energy EFT of QCD by saturat- 
ing the LECs, which has been termed resonance satura- 
tion [13-15]. We will discuss this issue briefly in Sect. El 
Here, we are interested in an extension of CHPT where 
these spin-1 fields are accounted for explicitly. While it is 
straightforward to construct the corresponding chiral ef- 
fective Lagrangian (as briefly reviewed also in Sect. El, the 
computation of loop diagrams is not. This is related to the 
appearance of a large mass scale in this EFT, namely the 
non-vanishing chiral limit mass of the spin-1 fields. This 
scale destroys the one-to-one mapping between the chiral 
expansion and the loop expansion, as discussed in more 
detail in Sect. 01 To be able to proceed in a systematic 
fashion, one has to be able to separate the contributions 
to loop diagrams originating from this large mass scale 
in a controllable and symmetry-preserving fashion. This 
problem also appears in the EFT when nucleons (baryons) 
are included (in fact, it has been analyzed first in this 
context [16]), and various solutions have been suggested, 
like heavy baryon CHPT [17, 18], subtraction schemes for 
the hard momenta [19], infrared regularization [20] or ex- 
tended on-mass shell regularization [21]. For the case of 
the vector mesons considered here an additional compli- 
cation arises - these particles can decay into Goldstone 
bosons and thus appear in loops without appearing in 
external lines. Here, we will present an extension of in- 
frared regularization that allows to treat such diagrams. 
In Sect. El we briefly review the intuitive approach of 
Ref. [19] how the contributions from the hard scale can 
be tamed. The more elegant infrared regularization [20] 
is introduced in Sect. In Sect. |BJ we discuss the new 
contributions to the Goldstone boson self-energy graph 
when the heavy particle only appears in the loop. The 
singularity structure of these loop graphs is analyzed in 
Sect. Based on that, we show how the infrared singular 
part can be obtained for such type of one-loop diagrams 
in Sect. El The method is then applied to the self-energy 
graph where the spin-1 field only appears inside the loop, 
see Sect. El The corresponding triangle graph is discussed 
in Sect. El Section ITT1 contains the analysis of the vector 
meson self-energy diagram with a pure Goldstone boson 
loop, which does not only involve soft momenta. As an 
application, we discuss the chiral extrapolation of lattice 
QCD data for the rho meson mass and related topics in 
Sect. El A brief summary is given in Sect. El Some tech- 
nicalities are relegated to the appendices. For other works 
on the problem of vector mesons in chiral EFT, we refer 
to [22-28]. 



2 Prelude: Vector mesons in trees 

In this section we show how vector mesons can be treated 
in chiral perturbation theory when only tree graphs are 
considered. This material is not new, but is needed to set 
up the formalism and to set the stage for the discussion 
of vector mesons in loops. The reader familiar with this 
material might skip this section. Also, our considerations 
are more general, they really refer to the coupling of heavy 



degrees of freedom to the Goldstone boson fields. Further- 
more, when talking about vector mesons, we really mean 
vector and axial- vector mesons (spin-1 fields). 

Our aim is to write down an effective Lagrangian con- 
taining the vector meson resonances explicitly. The word 
'explicitly' refers to the fact that the vector meson reso- 
nances are present implicitly in the Goldstone boson effec- 
tive Lagrangians through their contributions to (some of) 
the pertinent low-energy constants, see Refs. [13-15]. We 
want to state this on a more formal level. Assume that we 
have constructed a Lagrangian Cn es (R, U, v, a, s,p) where 
R are some resonance fields which might be the vector 
mesons for example. U collects the Goldstone bosons and 
v,a,s,p are vector, axial- vector, pseudoscalar and scalar 
sources. The latter also include the quark masses, s(x) = 
M, with M = diag(m M , mj, m s ) the quark mass matrix. 
The resonances are all very much heavier than the Gold- 
stone bosons (e.g. M p ~ 770 MeV), and therefore it is a 
consistent procedure for a low-energy effective theory to 
integrate out these heavy degrees of freedom by means of 
a path integration over R: 

JyiJljgiS d^xCRes(R,U,v,a,s,p) _ e iZ ind (U,v,a,s,p) ^ ^ 

By doing the path integral, a Goldstone boson theory 
(containing only U and the external fields) is recovered. 
Zind(U,v,a, s,p) may be called the generating functional 
induced by Cji es (R, U, v, a, s,p) through the path integra- 
tion. Such a step may be visualized in a Feynman graph 
by shrinking the lines symbolizing the i?-propagators to a 
point, and attaching contact terms (interactions between 
the remaining fields) to these points. 

What is the physical content of all this? Computing 
the i?-induced Goldstone boson field theory means that 
some interaction terms between the Goldstone (and exter- 
nal) fields are computed by the path integration over R 
and are therefore expressed through couplings of the res- 
onance to these fields and the resonance mass Mr, which 
can be measured in processes where the resonance occurs 
as an external state. The interactions generated in this 
way can be compared with the couplings determined by 
the LECs of the original Goldstone boson field theory. This 
will show us how important the resonance contributions 
to the processes under consideration are. We get a micro- 
scopic information on these processes which we could not 
achieve with the theory of the light fields alone. 

Numerically, however, one could as well work with the 
original effective field theory and simply include higher 
and higher orders. As should be clear by now, the reso- 
nance contributions are implicitly present in this theory, 
influencing the values of the LECs. This can be illustrated 
by considering a diagram with a resonance line through 
which a small (0(p)) momentum flows. The resonance 
propagator can be expanded in this case using 

w^Mi-w( 1+ m + -)' (2 ' 2) 

leading to an infinite series generating terms of an arbi- 
trary high order. Including the resonance field explicitly, 
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one takes care of all terms in this series and not just the 
first few terms of it. This can be advantageous. A nice ex- 
ample can be found in [29] , where the inclusion of vector 
mesons substantially improves the results for the nucleon 
form factors computed in that work. 

With this motivation, let us now concentrate on the 
vector mesons. Following [13], we will (first) use an anti- 
symmetric Lorentz tensor-field = —W V u. to describe 
the vector meson. This has six degrees of freedom, but we 
can dispose of three of them in a systematic way, for de- 
tails see Add. 1X1 and Ref. [13]. The spin-1 fields transform 
as any matter field under the non-linearly realized chiral 
symmetry, 

W MV (x) -» hW^xjh) , (2.3) 



where the compensator field h is defined via 

u(x) -> g R u{x)h) = hu{x)g\ 



(2.4) 



with gi is an element of SU(3)j, I = L,R and u 2 = U. 
The kinetic and mass term of the effective Lagrangian for 
vector mesons has the form 



Now we examine the interaction of the vector meson field 
with the other fields of the theory, especially the Goldstone 
boson fields. We actually already have interaction terms 
coming from the covariant derivative, but they are bilinear 
in Wfa,. For the simplest diagrams we want to consider, 
we need vertices with only one vector meson line attached 
to them, i.e. couplings linear in W^ v . Therefore we neglect 
the 'connection terms' in the following. From the philos- 
ophy of effective field theories, we are required to con- 
struct the most general interaction terms consistent with 
Lorentz invariance, chiral symmetry, parity, charge con- 
jugation and hermiticity. The building blocks with which 
this can be done are, in principle, at hand: the Goldstone 
boson fields (collected in U or its square root u), the co- 
variant derivative acting on U, the object x which 
contains the scalar and pseudoscalar fields (in particular, 
the quark mass matrix), the field W^u (and the pertinent 
covariant derivative V M acting on it) and the field strength 
tensor for the external fields and a M . For our pur- 
poses, it is more convenient to collect these blocks in the 
combinations 



£w = -iiVW^VpWn + -M*(W^Wn , (2.5) 



where 



W — W a T a 



(2.6) 



for an octet of vector mesons, and summation over the 
flavor index a = 1, ... ,8 is implied, and ( ) denotes the 
trace in flavor space. The pertinent covariant derivative is 



V fl R = d fl R+[r^,R] 



(2.7) 



It transforms as Wpv under the chiral group. Here, is 
the connection, 

= ^{v)[d^~i{v^ + d L )]u+u[d^-i{v^-d L )]v)) . (2.8) 

For the SU(3) case we consider here, the T a are the usual 
Gell-Mann-matrices which obey (T a T b ) = 2S ab as well as 
[T a ,T b ] = 2if abc T c , where f abc are the totally antisym- 
metric structure constants of SU(3). The mass My ap- 
pearing in Ea. Q2.5|) is strictly speaking the vector meson 
mass in the chiral limit. 

The vector meson octet we consider here may be given 
in matrix form: 



W 




K* 



(2.9) 



From the above Lagrangian, one can derive the propaga- 
tor, 



":r(w;„«fe))|o) = ^ 

I /r2 



d 4 k 



%k[x—y) 



(2tt) 4 M 2 - k 2 - it 



x - k ) + g^pkvka - g^Kkp - (p <-» v)] . 

(2.10) 



( p = iu'DpUu' = , 
a, = iw D^DvUv) , 



= „t^,t ± „^t 



X± = u'xu 
p± = uF L 



F± = uFLu^ ± v)F?;„u 



ux'u , 



This is better from a practical point of view because the 
so-defined objects all transform like under chiral trans- 
formations. This makes it easy to find chirally invariant 
expressions: Just take some of the above objects and put 
them inside a trace (...). This will then be invariant by 
the cyclicity property of the trace. Of course, one can fur- 
ther reduce these possibilities by imposing the other sym- 
metries mentioned above, and using the antisymmetry of 

What concerns power counting, it is clear that is a 
quantity of order 0(p) due to the covariant derivative D M 
giving one factor of momentum or external (axial-)vector 
source. Similarly, the other objects in the list are of 0(p 2 ). 

At order O(l), there is no term linear in W because 
W aa — 0. The lowest order interaction terms turn out to 
be of order 0(p 2 ) and read [13] 



r int 



iGy 
2y/2 



([ Ufl , Uv }wn- (2.11) 



This is more complicated than it looks, because both terms 
contain interactions with an arbitrary high (even) number 
of Goldstone boson fields. One must carefully expand the 
objects F^, and to obtain the vertex for a particular 
amplitude. It is now clear why vector meson singlets can 
be neglected: The sources to which they would be allowed 
to couple are (F^ v ) and ([u^, u u ]), but both traces are 
zero. 

Note that a discussion of the numerical values of Fy 
and Gy is given in Refs. [13,15]. Furthermore, the vector 
field formulation is summarized in App. El 
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3 Vector mesons in loops: Statement of the 
problem 

Difficulties arise when in a Feynman diagram lines rep- 
resenting a heavy matter field are part of a loop. The 
corresponding amplitude will in general not be of the chi- 
ral order expected from power counting. This has been 
noted many years ago when the nucleon was incorpo- 
rated in CHPT [16]. We already saw in the last chapter 
that power counting in CHPT is not as straightforward 
as in the purely Goldstone bosonic sector, and we already 
noted the reason for this, namely, the presence of a new 
large scale: The mass of the heavy matter field. In [16], it 
was shown that the parameters of the lowest order pion- 
nucleon Lagrangian already were (infinitely) renormalizcd 
by loop graphs in the chiral limit. There exists a mismatch 
between the loop expansion in h and the chiral expan- 
sion in small parameters of order 0(p). The loop graphs 
in general generate power counting violating terms, con- 
fusing the perturbative scheme suggested by CHPT. For 
example, a graph with dozens of loops might give a con- 
tribution as low as 0(p 2 ). Clearly, we will have to get rid 
of these power counting violating terms if we want to keep 
this scheme when including heavy matter fields. 

For graphs where the heavy field only shows up in in- 
ternal tree lines, the problem concerning power counting 
was not urgent, because the chiral expansion of the ampli- 
tude corresponding to such a graph at least started with 
the correct order. We saw this in the last chapter: We 
counted the vector meson propagator as O(l), and in the 
low-energy region, the momentum transfer variable t was 
much smaller than the square of the heavy mass, allowing 
an expansion of the propagator in the small dimensionless 
variable t/M v , which starts at 0(1). 



4 Soft and hard poles 

If the heavy particle line belongs to a loop, an integration 
over the four-momentum flowing through this line takes 
place. Due to the pole structure of the propagator, the 
integral will pick up a large contribution from the region 
where the line momentum squared k 2 w M v , with My 
the heavy mass (we keep the index V to remind us that 
we will concentrate on vector mesons, but the present dis- 
cussion is more general). This is the region of the 'hard 
poles' in the terminology of [19]. These contributions are 
of high-energy origin and clearly do not fit in a low-energy 
effective theory - they must be identified as generating the 
part of the loop integral which violates the power count- 
ing scheme, because this scheme would clearly be valid if 
these 'hard poles' were missing. The only hard-momentum 
effects for loop integrals in the Goldstone boson sector are 
the ultraviolet divergences, which are handled by dimen- 
sional regularization. 

Before discussing the infrared regularization method, 
for illustrative purposes we will shortly present the idea 
of [19]. We remarked that that the power counting vio- 
lating terms stem from the region k 2 w M 2 . Far off that 



region, for k 2 <C M 2 , it would be allowed to expand the 
propagator in the small variable k 2 /M 2 , as we did for 
tree lines in the last section. Of course, this is not allowed 
under a momentum integral which extends to arbitrary 
high momenta /c M . But it is too much of a temptation to 
do so, because in this way one destroys the 'hard poles' 
responsible for the power counting violating terms. Doing 
the expansion, treating the loop-momentum fc M over which 
one integrates as a small quantity, and interchanging in- 
tegration and summation of this expansion, one ends up 
with an expression which obeys power counting, because 
the hard poles are not present in any individual term of 
the series which one integrates term by term. Only the soft 
poles from the Goldstone boson (or perhaps also photon) 
propagators will be present in the individual integrals. 

Clearly, this is not the old integral any more, but a 
certain part of it, collecting only the contributions from 
the 'soft poles' - it is called the 'soft part' of the full inte- 
gral in the terminology of [19]. The remaining part, which 
was dropped by this procedure, collects the contributions 
from the 'hard poles' and stems from the high-energy- 
momentum region. It is argued that this part is expand- 
able in the small chiral parameters and, truncated at a 
sufficiently high order (depending on the order to which 
one calculates) it is a local polynomial in these small pa- 
rameters and can be taken care of by a renormalization of 
the parameters of the most general effective Lagrangian. 
The power counting violating terms are then hidden in 
the renormalization of the LECs, and the 'soft part', with 
subtracted residual high-energy divergences, is taken as 
the renormalizcd amplitude appearing as a part of the 
perturbation series. If this argument is correct, and the 
part of the full integral which is dropped is indeed ex- 
pandable in the small parameters, the 'soft part' contains 
all the terms of the full integral which are non-analytic in 
the expansion parameters, like terms of the typical form 
lnM«^//x, the so-called 'chiral log'-terms where M$ is the 
Goldstone boson mass, and \i is the renormalization scale 
(we will use dimensional regularization throughout). This 
becomes large if one lets the Goldstone boson mass go to 
zero. The physical interpretation is that in this limit the 
range of the interaction mediated by the Goldstone bosons 
becomes infinite, so that, for example, the scalar radius of 
the pion (or other hadrons) is infrared divergent in the 
chiral limit, which makes sense because it measures the 
spatial extension of the Goldstone boson cloud. 

We keep in mind that the 'soft poles' in the region 
where the loop momentum k <~ 0(p) are responsible for 
the terms non-analytic in the small parameters, like chiral 
log terms, and that all these terms obey a simple power 
counting. This follows from the argument given by [19]. 
If one ever encounters a power counting violating term 
non-analytic in the quark masses, for example, this would 
invalidate the above argument - such terms can not be 
hidden in the renormalization of the parameters of an ef- 
fective Lagrangian. 

We would like to demonstrate that arguments such as 
the one cited above are indeed more than just wishful 
thinking. The basic features of the argument can be ex- 
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hibited by a toy model loop integral (it is straightforward 
to generalize it - the formulas would only look a bit more 
complicated) . 

Consider the integral 



dk 



f(k) 



2m (k -a)(k- b) 



(4.1) 



where C is some contour in C \ {a, b} (enclosing the poles 
or not) and f(k) is an analytic function. Such an integral 
might e.g. be the fc°-integration over a full loop integral. 
The soft pole is called a and is of 0(p) while the hard pole 
is b. In fact, we will only need \a\ < \b\. In the spirit of the 
argument of [19], we destroy the hard pole in two steps 



dk f(k) 
q 2ixi (k — a) 

dk 



It 



1 
b 

f(k) ^ 



2m (k - 



k 



t 2 



+ 



soft 



Depending on the form of f(k), this will be ultraviolet 
divergent if the contour extends to infinity, but this diver- 
gence has nothing directly to do with the pole structure 
and can be handled by a regularization scheme. By the 
method of residues, / so ft is computed to be 



dk 



1 



(a — b) Jq 2m' (k — a) 



(4.2) 



We were allowed to sum the geometric series because \a\ < 
\b\. Therefore I so ft is just the first summand in the decom- 
position 



dk 



f(k) 



q 2iri (k — a){k — b) 
1 f f dk f(a) 



dk f(b) 



a — b \Jc 27ri (k — a) J c 2ixi (fc — 6) 

= -fsoft + -fhard , (4.3) 

which clearly separates the contributions from the soft and 
the hard pole, respectively. We have thus demonstrated 
that the expansion of the hard pole structures followed by 
an interchange of summation and integration indeed gives 
exactly the contribution from the soft pole. Note that the 
above argument can be iteratively used for multiple poles. 

A complementary approach would be to treat k as 
large, expand the soft pole structure in the small vari- 
able a and again interchange summation and integration, 
thereby isolating the hard pole contribution which is by 
construction analytic in the small parameter a. This may 
be called the 'regular' part of the loop integral. Again, 
divergences can show up which have nothing to do with 
the details of the pole structure (infrared divergences if 
the contour C encloses k = 0), but apart from this, the 
calculation analogous to the preceding one will give noth- 
ing but /hard- In this sense, the approaches of [19] and [21] 
can be said to be complementary to each other. 

We will see in the next section that the arguments 
of [19] can be refined or, as one should say, modified in a 



rather elegant way. In particular, the method described in 
this section does not always work because not all integrals 
converge in the low energy region, for more details on that 
issue, see e.g. [20]. 



5 Infrared regularization 

In order to find a more elegant way to separate the low- 
energy part of the loop integrals, we will examine a loop 
consisting of one 'heavy' propagator and one Goldstone 
propagator. We use dimensional regularization to handle 
the ultraviolet divergence of such an integral, i.e. we give 
a meaning for the notion of an integral in d dimensions, 
where d might be fractional, negative , etc. . 
We consider the scalar loop integral 



Iv<t>{l) 



d d k 



1 



(27r) d {{k-qf 



(5.1) 



Here q denotes an external momentum flowing into the 
loop, M<j, is the Goldstone boson mass and My is the 
mass of the heavy particle. The 'ie'-prescription, giving 
the masses a small negative imaginary part, should be 
understood here. We leave it out because it will play no 
role in the following discussion. 

There are two cases which have to be distinguished: 
1) the momentum q belongs to the heavy particle line 
of the loop, in the sense that this line is connected to 
external heavy particle lines (this would necessarily be the 
case if the heavy particle is a baryon, because of baryon 
number conservation). For the soft processes we consider, 
this would mean that 

q 2 - M 2 V = 0(p). 

The second case is that 2) the loop is connected only to 
Goldstone boson lines, which is the case for the vector 
meson contribution to the Goldstone boson self energy 
(Fig. m>) . This can not happen if the heavy particle line 
represents a baryon. Then 

q 2 = 0(p 2 ). 

We first investigate case 1). To this end we use a Feynman 
representation 



-4= / <l: 



1 



ab Jq (a(l — z) + bzf- 



to write 



Iv4> — 



d d k 
{2ttY 



dz 



l((k-qy-Ml){l-z) + {k*-M*)z] 



(5.2) 



Note from this expression that, for z = 0, the integrand is 
a pure 'soft pole', while for z — 1 it is a 'hard pole'. We 
see that the soft pole structure which we want to extract 
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Fig. 1. Self-energy graphs. Solid and dashed lines denote vec- 
tor mesons (heavy particles) and Goldstone bosons, respec- 
tively. While a) can be treated by the IR method of Ref. [20] , 
to deal with graphs of type b), the method developed here has 
to be used. 



is associated with the region z — > 0. After a shift k - 
k + q(l — z),the denominator of the integrand becomes 



where 



\e-A{z)f 



A(z) = M*C(z), 

C{z) = bz 2 - (b + a- l)z + a, 



(5.3) 



and the d— dimensional A:— integral can be done in a stan- 
dard way to give 



Iv<j> = ~- 



d-A 



(47T)i 



-r 2 



dz(C(z))i- 



(5.4) 



This will develop an infrared singularity as — > for 
negative enough dimension d. From the expression for 
C(z), we see that this singularity is located at z = 0, 
because in that region we will have 



\M V 



generating a part non-analytic in the Goldstone boson 
mass (remember that d can also be fractional). These find- 
ings are consistent with the above observation that the 
region z — > is to be associated with the soft pole struc- 
ture, coming from the Goldstone boson propagator. We 
come to the conclusion that to extract the soft pole con- 
tribution, we should extract the part of the loop integral 
proportional to nonintcger powers of the Goldstone boson 
mass (for nonintcger d). 

Bechcr and Leutwyler have proposed a way to achieve 
this [20]. They find that the decomposition of the loop 
integral into a part non- analytic in a (and therefore M^) 
and a part regular in a is given by 



I v<t> = I + R, 



where 



I = — 



d-A 



(4tt)3 

d-A 



-T 2 



R 



Mi 



(47T)i 



r 2- 



o 



'-dz 



2 dz. (5.5) 



From the above remarks it is clear that R, which contains 
the parameter integral starting at z = 1, will not produce 
infrared singular terms for any value of the dimension pa- 
rameter d. It will therefore be analytic in the Goldstone 
boson mass. On the other hand, the so-called 'infrared 
singular part' I is exactly the part proportional to nonin- 
teger powers of the Goldstone boson mass, for noninteger 
dimension d. Moreover, it is shown in [20] that this part 
of the loop integral fulfills power counting (as we would 
have by now expected for the part associated with the soft 
pole structure, see the discussion of the last section). We 
will not repeat the proof for these assertions here, because 
we are going to do a very similar calculation in the next 
section. 

The parameter integrals for / and R do not converge 
for d = 4. To give them an unambiguous value, a partial 
integration is performed to express them through conver- 
gent integrals and a part that is divergent for d = 4, but 
can be left away by an analytic continuation argument 
(analytic continuation from negative values of the param- 
eter d). For details, see Ref. [20]. This leads to an explicit 
representation of I and R for the case of four dimensions. 

The infrared regularization (IR) scheme can now be 
implemented by simply dropping the regular part R, ar- 
gueing that it can be taken care of by an appropriate 
renormalization of the most general effective Lagrangian. 
This can be done because it is analytic in the small pa- 
rameter (and in external momenta). The regular part 
contains the power counting violating terms originating 
from the 'hard pole' of the loop integral, which have now 
been abandoned with the dropping of R. We are left with 
the infrared singular part / which obeys power counting. It 
still contains a pole in (d — 4), which can be dealt with us- 
ing e.g. the (modified) minimal subtraction scheme. Hav- 
ing done this, we have achieved the goal of a finite loop 
correction where power counting allows to compute cor- 
rectly the order with which this correction will appear in 
the perturbation series. 

The method as presented here strongly relies on di- 
mensional regularization. In particular, the separation of 
the loop integrals into two parts with fractional versus 
integer powers of M$ for a fractional dimension parame- 
ter d allows the argument that chiral symmetry (or more 
specifically, the Ward identities) have to be obeyed by 
both parts separately. Therefore dropping one part of it 
(the regular part) is a chirally symmetric procedure be- 
cause it leaves us with a regularized amplitude that is 
again chirally symmetric for itself. It is important here 
that the scheme of dimensional regularization leaves chi- 
ral symmetry untouched, because the validity of the Ward 
identities docs not depend on the space-time dimension 
parameter. Of course one must deal with all loop inte- 
grals occurring in the perturbation series in the same way 
to keep the physical content of the theory unchanged. The 
regular part of any loop integral that one computes will 
have to be dropped. 
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If one lets d — > integer n, we will get terms proportional Here we have introduced the zeroes of C(z), 



to 



\M V \M V \ \M V 



In the complete expression, one will leave out terms of 
0(e), so that for an integer dimension, the expression for 
the infrared part / (up to 0(e)) may well contain a piece 
analytic in the Goldstone boson mass. Only after the sep- 
aration in the two pieces of different analyticity character 
has been done, it is allowed to let d approach an integer 
value. This is why dimensional regularization, permitting 
noninteger valued dimension parameters, is essential for 
this approach. Note that an elegant extension of infrared 
regularization to multi-loop graphs is given in Ref. [30]. 



6 Another case of IR regularization 

The last section, where the infrared regularization scheme 
has been introduced, dealt with the case that the momen- 
tum squared q 2 , coming from outside the loop, was of the 
same order as M v . This is the case, for example, for a 
diagram contributing to the self-energy of a nearly on- 
shell baryon with a Goldstone boson loop. This was the 
case considered in [20]. However, we have seen in the last 
chapter that vector mesons can occur as strongly virtual 
resonance states in processes where only soft Goldstone 
bosons or photons appear as external particles. If the same 
loop we considered in the last section is connected only to 
Goldstone boson lines, the momentum flowing into the 
loop will then be small. Since this is a Goldstone boson 
momentum, we must require that the 'regular part' which 
we want to drop from the regularized amplitude is analytic 
also in this parameter. 

Becher and Leutwyler introduce in their original paper 
the variable 



n = 



I 2 -Ml 



My 



2M <P M V 



(6.1) 



which is O(l) for the processes they examine, which corre- 
spond to the first of the two cases we distinguished in the 
last section. They consider the chiral expansion for fixed 
Q. But if 

\q 2 \ « My, q 2 = 0(p 2 ), 

the variable Q will be 0(jT x ) ! This already signals that 
this case will probably have to be treated differently. 

Let us first evaluate the integral directly, for d = 4 — e. 
We find (omitting terms of 0(e)) 

Iv * = -770 r f 2 - I) f dz ^ z x ^ z X2 ^~ 2 

(4tt)2 V 2 J Jo 



(47T) 

'-'A ■ ' . 1 



= 2A- 



16tt 2 16tt 2 
1 



dz ln(6(z — Xi)(z — x 2 )) 



16?r 2 



b + a-l , 
Xl < 2 = ^b~ ± 

and the standard notation 
1 



l) 2 -4ab 



4b 2 



(6.2) 



A 



My~ 4 

16tt 2 



d-4 - 2 ^)-l + D 



We have selected the mass of the heavy particle as a nat- 
ural choice for the renormalization scale \i here (this par- 
ticular choice also leads to suppression of higher order 
divergences that appear in the loop integrals and should 
thus be made). Furthermore we have used C(z = 1) = 
b(\ — x\)(\ — x 2 ) = 1 . To examine this further we write 
the expansions 



ab 



x\ 
1 

X2 



I- a (1 - af 
b b 2 



(6.3) 
(6.4) 



1-a (l-af 
showing that, for b — > 0, x 2 behaves like 

b + a-l 

Remember that in the case we consider now, a and b are 
both small variables from their definition, of chiral order 
0(p 2 ). Using these expansions, and the relation 

C(z = 0) = a = bxix 2 , 

we see that Iv<j> contains a term non-analytic in the Gold- 
stone boson mass 



1 



iv<p 



16tt 2 



xi ln(a) + . . . , 



and that it is analytic in the second small variable b (for 

We want to check the power counting for this case: The 
non-analytic terms are proportional to xi, whose expan- 
sion in a and b starts at order 0(p 2 ). This is the expected 
chiral order for the integral: The loop integration in 4 
dimensions gives 4 powers of small momentum, whereas 
the Goldstone boson propagator gives —2. The hard pole 
structure is of order 0(1) here, because the vector meson 
appears just as an internal resonance line. So the power 
counting for the non-analytic terms is fine, as expected. 

We remark that these non-analytic terms are also pro- 
duced when one proceeds after the prescription of [19]: 
Expanding the hard pole structure and integrating term 
by term, one gets x\ as the coefficient of the lna-terms, 
order by order. 

Now we want to do infrared regularization. But now 
we remark an important point: We can not simply take 
over the formulas of the last section. We note that the 
expression for / and R will contain pieces non-analytic 
in the other small variable b, which is not small in the 
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case treated in the last section. But the original integral 
does not contain such a non-analyticity in b. We conclude 
that the extension of the parameter integrals to z > 1 
is responsible for this non-analyticity: Somewhere in that 
region, we must catch up a pole for b — > 0. So we have 
the problem that, loosely speaking, the 'regular part' is 
not regular. We will have to modify the method of [20] 
somehow, and find a way to separate off the terms non- 
analytic in the small variable b. To do this, we will have to 
examine the nature of the singularity we encounter here. 
This will be done in the next section. 



7 Singularities in parameter space 

We will consider an integral which has exactly the same 
features as the one we need, but is a bit simpler. Let 



Jo 



dz (z + a)~ 3 (l + bzY 



(7.1) 



where a and b are again small parameters in the sense 
that \a\, \b\ <§C 1. This will develop a singularity at z = if 
a — > 0. What if we extend the integration to infinity? To 
examine this we first change variables, 



1 

z= -, 
u 



and compute 



u du (1 + au) 2 (u + b)' 



(7.2) 



This shows the close similarity between the 'a-singularity' 
at z — and the '6-singularity' at u = 0. If we extend 
the integration to z —* oo, we will pick up a non- analytic 
contribution from u = 0. A 'regular part' defined as 



R = 



dz (z + a) ^(l + bz)~ 



(7.3) 



will be analytic in a, but not in b. This is the situation 
encountered in the last section, for a slightly different in- 
tegral. How do we get rid of this non-analyticity? 

Remembering what we have learned so far, we know 
how we can get rid of certain non-analytic terms: We must 
'destroy the poles' by expanding the pole structures and 
integrating term by term. Let us do this: 



dz(z + a)-?(l-^bz + ^-(bz) 2 ±... 

. 3 ^ r(-i) (bz) m 

v ' ^T-i-m to! 



I 



m=0 



~ n-D b™ r 1 



z + a) 2 z". 



Here we were allowed to interchange summation and in- 
tegration without changing the value of the integral, be- 
cause in the interval of integration the expansion of the 



integrand is absolutely convergent when a,b are smaller 
than 1. 

The dependence on a now resides in the simpler inte- 
grals of the coefficients in the expansion. The idea is now 
to find the 'infrared singular part' of each of these coef- 
ficients. If the sum of these infrared singular parts con- 
verges, it must be the infrared singular part of the full 
integral /, because the expansion of the pole structure 
that we have performed did not change the integral. 

To find the infrared singular part of 



= / dz(z + a] 
Jo 



d „m 



(7.4) 



with an arbitrary parameter d (possibly negative) and an 
integer to, we follow the method of [20] and scale z = ay 
to find 



d+m+1 



dy (1 + yfy 



(7.5) 



We would now like to take the upper limit of the integra- 
tion to infinity. Whenever the extended integral converges, 
it gives 



Jo 



dy (1 + y) d y 



d ni m 



r(m+l)r(-d- (m + l)) 



There will be a convergence problem if to is large enough. 
But the divergence comes from large values of z and has 
nothing to do with the infrared singular part. To see this 
in detail, let K ^> 1 , m > 0, and do a partial integration: 



[ K dz (z + a) d z m - K m ( K + a) d+1 
Jo a + 1 



-f 

Jo 



d- 
m 



d+l 



(z + a) 



d+x z m - x dz . 



Doing this m times, and dropping all terms proportional 
to (K + a) d+n with n E N , because these are clearly 
expandable around a — and will therefore not contribute 
to the piece non-analytic in a, we end up with the same 
result as above, namely, the 'infrared singular part' of I m 
is always 

r _ n d+ m +i r(m + IJTj-d - (to + 1)) 

I m ,m - a . (7.6) 

This form clearly shows the character of the infrared sin- 
gular part as being proportional to noninteger powers of 
a for noninteger parameter d. 

We have now performed the relevant step to find the 
infrared singular part of each coefficient in the expansion 
of /. We insert this in the series, setting d = — | for defi- 
nitencss: 



E 



n-h) 



dz (z + a) ?z m 



E 



T(-i) r{m+l)r{\-m)(ab) m a-? 



-o r (- 



to) 



TO! 
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One can easily sum the series, 

T 4 V"^ 



'IR 



v m— 

4 / 1 



(ab) m [ (m + 1)- 



1 



/E\(l-ab) 2 2(1 -ab) 
2 1 + ab 



^/E{l-ab) 2 ' 

The reason why we have selected the value d = — § is that 
it is very easy to compute the integral I directly. The part 
non-analytic at a — > can be read off and is the same 
as the result just given. We have checked this for various 
other values of the parameter d. 

In a 'naive' application of the IR method, one would 
be tempted to simply take 



BL 



dz (z + a) 3 (1 + bz) 2 



(1 - ab) 2 



1 + ab 



-2Vb 



which contains a part non-analytic in the second small 
variable b. This is the part we have separated off by our 
procedure. 

We note that this '6-singular' part can be extracted by 
proceeding in exact analogy to the steps just performed: 
Expand the other pole structure, proportional to a power 
of (z + a), in the integral R, Eq. (|7.3|) . The expansion is in 
powers of a/z, which is smaller than one in the pertinent 
interval of integration. 

We can now give the result for arbitrary d. The general 
'a-singular part' is 



/IR = \ r{-d)r{d + i-m) {ab) a 

m— / 



m d-\-l 



(7.7) 



while the '6-singular part' is 



j h= J- r (d+l)r(m-2d-l)(abY (? 



m=0 



n-d) 



ml 



In the case d = — I , the last expression yields indeed 



\ " 



m=0 v 



(m+ l)(ab) m Vb 



-4Vb 



{1-abf 



which is confirmed by the result of the direct calculation. 
We have 

/ B L=Im + /b, (7.9) 

where the second part is the one we do not want in a 
genuine regular part (it will appear in R because the non- 
analytic behaviour for b — > is not present in the original 
integral /, with which we started). 

What we will need in the next section is the result for 
d — — e, where, as always, e is considered as sufficiently 



small to allow for the neglecting of terms 0(e 2 ). In this 

case, 



E 



T(l - e )r(-(m + l) + £ ) 
r(e)f(l - m - e) 



{ab) m a 1 - e . (7.10) 



After some /^-function algebra, one gets for the sum 
^— ' mm + 1 

771—1 V ' 

= a(-l-e + e H a))-ea±a-^-)(aby 



+0(e 2 ). 



The series can easily be summed 



In(a)+e( ^—^ ) ln(l-ab). (7.11) 



Please note that the last term is expandable in b, and that 
we have left out terms of 0(e 2 ). 

With the very same method, we can also compute the 
'6-singular part' for the case d = —e. The result is 

r.\ ) (]n(l - oft) - 1). (7.12) 



We have checked that 7bl, calculated with the method of 
Becher and Leutwyler, is again the sum of the 'a-singular 
part' and the '6-singular part', like it was the case for 
d = -|. 

We have now achieved a method which allows to split 
integrals of the form of 7bl into an 'infrared singular part', 
behaving non-analytically as a — > 0, and a part showing 
such behaviour for 6^0. Note that this decomposition is 
unique, and that both parts are of a different analyticity 
character (for fractional d) , concerning the small parame- 
ters a and 6, respectively. In the next section, we will sec 
how this method can be applied to the scalar loop integral 
Iv4>- 



8 Corrected infrared singular part 

The integral we need for the calculation of Iy^ is 
dz (b(z — xi)(z — x 2 ))^ 2 - 

Extracting a factor 

(_ te2 )#- 2 = (l-(a + &) + ...)*- a , 

from the integral, which is expandable in a and 6, the 
remainder is of the form of / treated in the last section, 
because x\ and small parameters of 0(p 2 ) : 

I={-bx 2 )i- 2 I' 
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where 

V = ( dz (z + (- Xl ))i- 2 (1 + {-x 2 )- l z)i- 2 . 
Jo 

Doing the appropriate substitutions in Ea. lj7.12jl . the in- 
frared singular part of I becomes 

ilR = x 1 +ex 1 -^x 1 \n(a)-^(x 1 -x 2 )\n^l-^J, (8.1) 

where we have used a — bx\x 2 . For completeness, we also 
rewrite the '6-singular part': 

4 = ~{x 2 -xi)Jn(6x|) 



(8.2) 



As a check, we add it to the infrared singular part: 
Im+h = z (l+e-^ln(a)J-^(x 1 - X2 )ln(^y (8.3) 

Here we have used the notation 

x\ + x 2 



as in Ref. [20]. Again, the sum of the 'a-singular part' 
and the '6-singular part' is the result for the integral 



dz(b(z — x\){z — x 2 )Y 



when computing it utilizing standard IR. But the correct 
infrared singular part for our case is only a certain part of 
it, namely, Zir,. 

The part which must be split off here, lb, vanishes if 

xi = x 2 => (b + a - l) 2 - 4ab = 0, 

(see Ea. (j6.2jl ). which is the case for 

q 2 = (My ± M ) 2 EE q 2 ± . 

We cannot trust our procedure for q 2 > My. Therefore 
the value q 2 _ should be considered as the point where the 
standard infrared singular part of [20] and the represen- 
tation given here, i.e. Tir, can be joined together. 

We emphasize that the kind of argument we have given 
here is completely in the spirit of the method of Becher 
and Leutwyler, in that we examined the analyticity prop- 
erties of the parameter integrals for a general dimension 
parameter d. 

The calculation of the infrared singular part of the 
scalar loop integral Iv<p can now be completed: 



d-A 



rlR 
1 V(p 



= 2xiA 



1 



16tt 2 



x\ — x\ ln(a) 



x \ 

-(xi - x 2 ) mil- — 
x 2 



(8.4) 



Terms of O(e) have been omitted. We claim that the dif- 
ference 

R' = Iv(j> — Iv\ 

is the appropriate regular part. This means that it is ex- 
pandable in the small parameters a and b around zero. 
The proof consists of two observations: 

1) Both Ivij, and Iy^ contain the same terms non- 
analytic in a, namely, 



1 



167T 



: xi ln(a). 



The difference is therefore expandable in a. 

2) Ivtp was expandable in b from the start, whereas Iyi 
is expandable in b by construction. Therefore the difference 
is of course also expandable in b. 

Moreover, R' is unique, because we extracted exactly 
the part of Ivtf> proportional to fractional powers of a for 
fractional dimension parameter d. We conclude that R! is 
a well-defined regular part, and that it can be absorbed in 
a renormalization of the LECs of the effective Lagrangian. 

We add the remark that the expansion of Eq. lj8.4j l is 
reproduced by using the procedure of [19], i.e. expanding 
the 'hard pole structure' and interchanging summation 
and integration. We have checked this to order 0(p 8 ), but 
a formal proof that it will give the same result to all orders 
is still missing. It seems that both procedures are indeed 
consistent (remember, however, the remarks made at the 
end of Sect.0J). This means that the 'low-energy-portion' 
of loop integrals is, in this sense, unambiguous. This result 
is not really surprising: From the arguments of Sect.0] it 
is seen that an integral like 7 so f t (see ea. (14.21) 1 is a pure 
'soft pole' integral, i.e. only involving the pole structure 
associated with the Goldstone boson propagator, and thus 
having no regular part, while an integral like /hard is reg- 
ular in the Goldstone boson mass, and does not contain 
fractional powers of M<p for any choice of the dimension 
parameter d. 



9 Goldstone boson self-energy 

We are now ready to compute the vector meson contri- 
bution of the Goldstone boson self-energy. We consider 
first the novel type of diagrams where the heavy mass line 
appears in the loop, see Fig.^i. We get for this amplitude 



12iGy b 



(9.1) 



using the notation 
j , , ' ddk 



q 2 k 2 - {q ■ kf 



{2iT) d (k 2 - M 2 + ie)((k - q) 2 - M 2 + ie) 

This integral may be decomposed in a linear combination 
of scalar loop integrals by standard techniques: 



Is — Ccf>Icf> + Cyly + Cy^Iytf), 



(9.3) 
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where the coefficients c, are given by 
q 2 + M 2 



c$ = 



Ml 



c v 



Ml 



Ml 



(q 2 -M 2 + M 2 ) 2 



CV<p 

while the scalar loop integrals are defined as 

Ml 



Iv 



/d k i 
J2^Y k 2 -M v +ie 



2M V \ 



(9.4) 

■ln(a), (9.5) 
(9.6) 



and the scalar loop integral Iy$ is defined in Ea. (|5.1|l . We 
repeat the remark that we use \i = My for the renormal- 
ization scale. 

What power would we like to have for this self-energy 
amplitude ? We have one loop integration, two vertices 
of order 0(p 2 ) and one Goldstone boson propagator. The 
vector meson propagator is counted as O(l) here, since the 
vector meson line is not connected to any external heavy 
particle lines. So we end up with the 'expected' power 

D(S) = 4 + 2x2 - (2 + 0) = 6, 

using Eq. (|1.2(l . The word 'expected' was used from a naive 
point of view, because we are already sophisticated enough 
to expect that the 'hard pole structure' associated with 
the vector meson propagator will give the loop integral a 
high energy contribution that spoils the power counting. 

Indeed, using the decomposition in scalar loop inte- 
grals and the expansions of the variables x\ and x-i given 
in eq.(3.10) and (3.11), respectively, it is straightforward 
to see that Is contains the following terms which violate 
the power counting: 



Is = \m v [ \(6b-2b 2 +6aby 



1 



16tt 2 



ab 5,-, 

b 

2 6 



where the dots stand for terms satisfying the power count- 
ing, i.e. they are of order 0(p 6 ) or higher. 

To find the 'infrared singular part' of Is, it is sufficient 
to find the infrared singular part of each of the scalar loop 
integrals. This is because the coefficients c; do not contain 
any fractional powers of for any dimension parameter 
d. 

The infrared singular part of Iy^ has been computed 
in the last section. The integrand of Iy is a pure hard pole 
structure without any dependence on a small parameter 
like a or b, and will therefore not have an infrared singular 
part. Finally, the integral 1^ is proportional to a fractional 
power of M^, as a direct calculation using dimensional 
regularization shows, so it has no regular part (it does not 
contain a hard pole structure which could be expanded). 
The 'infrared regularized' self-energy amplitude is thus 



1 E 



YliGy bf 



cv-ply^) 



(9.7) 



where Iy\ is given in eq. 18.4J1 . Using the expansions of the 
Xi and ln(l — y) = —y — y 2 /2 — . . . (\y\ < l), it can easily 
be checked that the 'infrared regularized' amplitude obeys 
power counting. 

Before we go on and apply our modified version of 
infrared regularization to other graphs, we want to men- 
tion one more thing. Using the vector field approach, cf. 
App. [5] the self-energy graph leads to the expression 



YliGy gabji 



M 2 F 4 



£ ■ 



where now 



I' r 



d d k 



k 2 (k 2 q 2 -(k-q) 2 ) 



(2ir) d (k 2 



M 2 +ie)((k~q) 2 -M 2 + ie)' 



Subtracting the amplitude computed in the vector field 
approach from the amplitude computed in the tensor field 
approach, we get 



A(W) - A(V) = 



IZiGy cab- 

MyF 4 



d d k 



{M 2 -k 2 ){k 2 q 2 ~{k-q) 2 ) 



(2Tr) d (k 2 



M 2 F A 



M 2 
d d k 



■ie){{k 
k 2 q 2 - 



-qf-M 2 
(k-q) 2 



■ie) 



(27r) d (k - q) 2 ~ M 2 + ie 



But this is the same result as one would get for a self- 
energy diagram where the vector meson line is replaced 
by a contact term interaction 



8M 2 



leading to a four-</> interaction 



G 



M 2 F A 



This confirms the 'duality' between the vector and the 
tensor field approach [15]. Note that this result does not 
depend on the regularization scheme - it was derived with- 
out even computing the integrals. Since the difference of 
the amplitudes is a pure 'soft pole' Goldstone boson loop 
diagram, we see that the 'hard part' (i.e. the regular part) 
is representation independent. In particular, both descrip- 
tions produce the same power counting violating terms. 



10 Triangle graph 

There is one more one- loop diagram of 0(p 6 ) where the 
vector meson line shows up as a loop line: the triangle 
diagram of Fig. [2J One gets the following expression for 
the triangle graph: 



Aa = 



GG 2 y 

F 4 



ab3 



[ fab 

V3 J 



i A> 



(10.1) 
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Fig. 2. Triangle graph as it contributes e.g. to the pion vec- 
tor form factor. Solid, dashed and wiggly lines denote vec- 
tor mesons (heavy particles), Goldstone bosons and external 
sources (fields), respectively. 



where the integral is 

f d d q 

-ir A {p,p + k)=J j^-y 

((2p + k) - 2q) T (p ■ (p + k)q 2 - (p ■ q)((p + k) ■ q)) 



(q 2 -M 2 )((q-p) 2 ~M 2 )((q-p-k) 



Ml 



(10.2) 

A decomposition of 1\ as a linear combination of scalar 
loop integrals is given in App.lDlffor the Goldstone boson 
momenta on mass shell). What concerns us here is the 
question how the infrared singular part of this integral can 
be obtained. The decomposition in scalar loop integrals 
contains an integral we have not yet treated, namely 



Iv<t>4>(p,P + k) 



d d q 
(27r) d 



(«* - M 2 )((q p) 2 M 2 )((q - (p + k)) 2 - M|) ■ 

(10.3) 

Fortunately, it is possible to reduce the problem of finding 
the infrared singular part of this integral to the case we 
have already examined. The procedure can in full gener- 
ality be found in section 6.1 of [20]. We show how this 
works in the above example: Introducing one more Feyn- 
man parametrization, we write Iy </>(/> as 



d d q 



(27r) d q 2 - M 2 



v Jo 



dw x 



d 



dM 2 (1 - w)((q - p) 2 - M 2 ) + w((q - (p + k)) 2 - MT) 



d d q 



(2tt)<* q 2 - M 2 Jo 
d 1 



dw x 



dM 2 (q-(p + wk)) 2 - (M 2 - k 2 w(l - w)) ' 

The momentum integral is now of the form of Ivcj>i with 
the operator 



Ai..,= l dw "(...) 



(10.4) 



acting on it. We can insert our result for J^, with the 
substitutions 

M 2 M 2 -k 2 w{l~w) _ , 



M 2 



M 2 



a , (10.5) 



b = 



(p + wk) 2 



My 



M 2 



b' 



(10.6) 



Please note that the external Goldstone boson momentum 
is now called p instead of q. Also note that the new vari- 
ables a' and b 1 are also of 0(p 2 ), which allows to take over 
the treatment of infrared regularization presented in the 
foregoing sections. 

We must show that the operator A does not disturb 
the properties of infrared singularity and power counting. 
The 'dangerous' part of this operator is the derivative with 
respect to M 2 , since it changes the chiral order. It is clear 
from the above definitions that (for fixed k 2 ) 

d d_ _ d_ 
dM 2 ~ da~ ~ da/ ' 



Mi 



We know from the derivation of the infrared singular part 
that it can be written in the general form 



(a')*" 1 EE' 



m— n— 



with some numerical coefficients c mn that depend only 
on the dimension d. For d — > 4, this gives the correct 
order 0{p 2 ) for Jy^. Letting the operator A act on this 
expression, we get 

I? M = AI?Ja',bf) 



dw 



J d 

M 2 da 



;((«') 



oo oo 

' i7 'EE'- "'" ,// " 

771—0 n— 



' m=0n=Q v 



+ m | c mn a' m b' n 



(10.7) 



This shows that the expansion of the thus defined infrared 
singular part of Iv<p<p starts with M d ~ 4 , as one expects 
for such an integral by simple power counting. We learn 
from the last expression that, in principle, it is sufficient 
to know the chiral expansion of Iy^ to arrive at the chiral 

expansion of /K,.. The only problem for practical calcula- 
tions is that the parameter integrals over w are not at all 
of a simple form, because a',b' and therefore also x\ and 
x' 2 , defined in analogy to Eq.(3.9), are nontrivial functions 
of w. 

For d -» 4, I 1 ^ is 



- (xi -4)ln(l-|)j) = 

ti^(2yi\-T^(yi-yiH*') 

-^-(w-».)Mi-i)- ja2 ^ > ) > ), 



where we defined 
d_ 



d , 1 , .. 



b' -a'-l 



y/{V + a' - l) 2 -Aa'b' , 
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The singularity structure of Iv<j><j> is richer than the one 
of Iv(f> because by the definition of a', a term like ln(a') 
not only contains the infrared singularity for — > 0, but 
also a cut for k 2 = t > 4M?, which is associated with the 
two-Goldstone boson production threshold. 

The important point is that, having the prescription 
for Ivrtn we can find the infrared singular part of any loop 
integral where a small momentum of order 0(p) flows 
through the heavy particle line(s) (this is the case we 
have treated in this paper) or where a nearly on-mass- 
shell heavy particle is involved (in which case the results 
of [20] can be used directly). The principle is now well- 
understood, but practical calculations will be difficult for 
complicated diagrams, because one needs a parameter in- 
tegration for every pair of propagators which are combined 
to one (parameter-dependent) pole structure. An exam- 
ple for this has been shown in the treatment of Ivtt><t>- As 
another point, the decomposition of a complicated loop 
integral involving a lot of vertices in a decomposition in 
scalar loop integrals will be lengthy and complicated. But 
these are no conceptual problems any more. The concep- 
tual problem of the power-counting violating terms has 
been solved by dropping the regular parts of all loop in- 
tegrals, retaining the infrared singular parts that stem 
from the region where the loop momentum is O(p). As a 
further consequence, many diagrams, namely those where 
the loops are formed of heavy particle lines only, can be 
dropped from the start because the respective loop inte- 
grals will only contain 'hard pole' structures and do not 
lead to an infrared singular part. Using this scheme, we 
can proceed and treat all kinds of diagrams where heavy 
resonances (not only vector mesons) occur in loops, and 
whose momenta are to be counted as either nearly on-shell 
or 0{p) by the perturbative scheme of power counting. We 
finally remark that the treatment of the vector meson loop 
graphs in the analysis of the nucleon electromagnetic form 
factors performed in [29] is consistent with the procedure 
we have established. 



11 Vector meson self-energy graph 

In the last section we have considered diagrams where the 
vector meson shows up as a strongly virtual intermediate 
state, with a small momentum flowing through the vector 
meson line. We found the 'infrared singular' part of the 
corresponding amplitude, and we saw that it was neces- 
sary to modify the method of [20], where all particles in 
the intermediate state where considered as being close to 
their respective mass shell. It is now natural to ask: What 
happens if the 'light' particles (the Goldstone bosons) are 
far from their mass shell ? In principle, the 'hard mo- 
mentum part' of the pionic intermediate states has been 
integrated out in the effective theory. But in analogy to 
the case of the treatment of the 'heavy' vector meson in 
the last chapter, it might be useful to take these degrees 
of freedom into account in a systematic fashion, because 
in this way one sums up (infinitely many) higher order 
graphs. We will encounter such a situation in the follow- 
ing section. 



11.1 One more case of IR regularization 

In Fig. |21 we show a graph contributing to the vector me- 
son self energy. In the case where there is a 'small' (0(p)) 
momentum flowing through the vector meson line, the 
corresponding amplitude would be a homogeneous func- 
tion of small parameters (external momentum and quark 
masses) , since the large scale (in this case, the vector me- 
son mass) does not show up in the loop line propagators 
and thus cannot produce a 'hard pole' contribution. The 
loop integral is the same as in the Goldstone boson sec- 
tor, and therefore has no 'regular part'. When computing 

P-k 

P 



k 

Fig. 3. The vector meson self-energy diagram with a pure 
Goldstone boson loop. Solid (dashed) lines denote vector 
mesons (Goldstone bosons). 

self-energy contributions, one is usually interested in the 
case where the external momentum P is close to the mass 
shell of the corresponding particle. This leads to the ap- 
pearance of the large scale in the denominator of the inte- 
grand of the loop integral, and we expect a power-counting 
violating contribution stemming from a 'hard pole' of the 
integrand. The integral will develop a regular part in the 
terminology of Bechcr and Leutwyler. 

We start our analysis for the case P 2 S> Mj with the 
question: What is the 'soft part' of a diagram like Fig. 4.1? 
That is, which region of the loop momentum integration 
produces the infrared singular part? Obviously, the case 
where both Goldstone boson propagators are of 0(p~ 2 ) 
is excluded by four-momentum conservation at the two 
vertices. The region where both Goldstone boson lines are 
far from their mass shell is a pure 'hard- momentum effect' 
and thus belongs to the 'regular part'. The soft part can 
only come from the region of the loop integration where 
one line is soft (i.e. it carries an O(p)-momentum), and the 
other Goldstone boson line carries the large momentum P 
and is thus far from its mass shell. 

As an illustration of this argument, let us first try to 
extract the 'soft pole contribution' following the method 
of [19]. We examine the scalar loop integral 



J (2nr(k 2 -M 2 )((k~P) 2 -M 2 y v ' > 

Following the above reflections and the receipt of [19], we 
treat the momentum of one line as 'soft' and expand the 
propagator associated with the other line. Then we inter- 
change integration and summation of the series, thereby 
'destroying the hard pole': 

. r d d k l i ^ (2P ■ k + m 2 - k 2 ) n 

(2ir) d k 2 - M 2 P 2 ^ {P 2 ) n 
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d d k 



(2ir) d k 2 -M 2 P 2 ^ (P 2 ) n 



oo 
p2 X/ 



{2P-k) n 



L^i (p2)n+l 
n=0 y ' 

TT-^soft- 



d d k {2P-k) n 

{2-K) d k 2 - Ml 



The factor of \ appears due to the fact that the full soft 
part also includes the part where the other line (with mo- 
mentum P — k) is considered as 'soft'. Of course this part 
is equal to the above result due to the symmetry of the 
graph. We note further that it is not legitimate to resum 
the series in the above result to get 



I' = i 



d d k 



1 



(2ir) d (fc 2 - Mj)(P 2 - 2k ■ P) 
d d k ro ° 



(2ir) d 



dz 



1 



[(k 2 ~ M 2 )(l - z) + z((k - P) 2 - M 2 )} 2 
d d k 



(2ir) 



dz 



1 



[{V-Ml){l-z)+z{{k-P?-Ml)r 

This contains a 'hard pole' contribution and will not sat- 
isfy the power counting scheme, which requires the scalar 
loop integral to be 0(p d ~ 2 ), because only one Goldstone 
boson propagator is booked as 0(p~ 2 ), while the other 
Goldstone boson must be far off its mass shell (its mo- 
mentum must be of the order of P by momentum con- 
servation). To repeat, this power counting is strictly valid 
only for the 'soft pole' part of the integral, which we have 
identified as I so ft- 

The alert reader will note that the result 7 so ft corre- 
sponds to a series of tadpole graphs, involving only one 
Goldstone boson propagator. This can of course not be 
the whole story, because the amplitude of Fig. [3] has an 
imaginary part due to the production of two Goldstone 
bosons in the intermediate state, while the tadpole sum 
does not have such an imaginary part. In order to take 
only / so ft as the regularized amplitude, one would have 
to write complex coefficients in the effective Lagrangian, 
which we do not want. A direct calculation of the full 
scalar loop integral shows that the imaginary part does 
not satisfy the power counting mentioned above. This is 
related to the fact that for large P 2 of the heavy external 
particle, the Goldstone bosons produced in the decay of 
this particle are not to be considered as 'soft'. Below the 
threshold, we have P 2 < 4Mj, so P 2 can not be consid- 
ered as being very large compared to the scale Mj in that 
region, and we would have to take the full integral I as the 
soft part, and not / so ft- This phenomenon of the 'missing 
imaginary part' is consistent with the findings of Ref. [25], 
where this was noted using the Heavy Vector Meson ap- 
proach. We will not discuss this further at this point and 



turn to the scheme of infrared regularization. Doing the 
usual steps, we obtain 



I = - 



£^l(P 2 )i- 2 I dz(D(z))l- 2 



where 



D(z) = z 2 -z + 



P 2 



(11.2) 



(11.3) 



Motivated by the remarks made in the last paragraph, we 
will consider the case that P 2 > 4M|. This is fulfilled 
for the case we are interested in, where P 2 is close to the 
physical vector meson mass squared, and is the mass 
of the particles we consider as Goldstone bosons. 

Obviously, fractional powers of are produced in the 
parameter regions where 

z 2 -z = 0^z = 0orz=l, 

corresponding to the fact that either one or the other 
Goldstone boson line in the loop carries soft momentum. 
In accord with the procedure of Sect. EI (cf. Eq. H6.2I) ). we 
introduce the zeroes of D(z), 

dx,2 = ^( 1= F^), 



<T= \ 1- 



AM 



4> 



p2 



(11.4) 



Note that a € R and 

< di - d\ = a < 1. 

We can simplify our analysis by 'folding' the parameter 
interval symmetrically, 



dz(D{z)) i - 



2 I dz{D{z))i 
lo 



#-2 



allowing us to expand the pole due to the zero c?2 > h 



dz(D(z)Y- 



dz(z — di)z 2 (z-d 2 )- 



= 2 



771=0 



{~d 2 )i- 2 -™ 1 n ~ 2 ~ l) 



ml r(| - 1 - rn) 



x / (z-di) 
Jo 



i~ 2 z m dz. 



We did not yet change the value of the parameter inte- 
gral. To find the 'infrared singular' part of the parameter 
integral in the last line, we note that d\ is proportional 
to -M| and of 0(p 2 ), and shift the integration variable to 
write 



(z- d 1 )2- z z m dz = z*- 2 (z + d x ) m dz 

(i J-di 

z*- 2 (z + d x ) m dz. 
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Terms proportional to dj 2 will only be produced by the 
first term on the right-hand side (remember m 6 N). 
Scaling the variable of integration with d±, it takes the 
form 



z^ 2 (z + di) m dz 



(-lr+^-dx)- 



-1+m 



ti- 2 (l-t) m dt 



= ( -iru)«-^ r( '-f (m+1) . 

v ; 1 ; T(f + m) 

where we substituted z = —td\. 

The 'infrared singular part' of I is thus 



E 



2r(2- f)(p 2 )f- 2 

(4tt)* 

(-i) TO (di)*- 1+w, (d2)i 



(11.5) 



i (nf-i)) 2 



r(f - 1 - m)r(f + m) 



This expansion starts with d-f ~ M? and obeys low- 
energy power counting. The series could be summed up, 
but this is not necessary. Reviewing what we have done so 
far, it becomes clear that we have just selected a certain 
range of integration which produces the fractional powers 
of M<f,. This step may be symbolized by 



/= / dz{...) 
Jo 



di rl 

dz{. ..)+ dz(. . .) 



dz(. ..)=/] 



IR- 



(11.6) 



Applying this to I with d = 4 — e, we find 



I m = 4diA- 



16tt 2 



., , -2di+ln(a)+oTn( )-2aln(a] 



1 - (7 



while the 'regular part' is 
1 



(11.7) 



I-h 



IB 



(2-4di)A- 



16tH 



-(l-2di)+2crln(o-)-i7r£7 ), 



(11.8) 

which is indeed expandable in Mj for P 2 > 4M^. We 
have again used My for the renormalization scale, and 
the variable a defined in Eq. Q5.3|l . 

It may be checked by expanding d\ and cr in powers 
of Mj that the infrared singular part indeed satisfies the 
power counting rules, and also that 

Im. = ^soft- 

We have already remarked in the last chapter that the low- 
energy part of a loop integral is unambiguously defined in 
this sense. 

The imaginary part of the scalar loop integral I is 

— icr 

167? 



whose chiral expansion starts 0(1) and therefore does not 
obey the power counting rules. But it cannot be sub- 
tracted from the full amplitude since it is not real. The 
corresponding width of the vector meson due to its possi- 
ble decay into a pair of Goldstone bosons cannot simply 
be neglected. In principle, one should give the denomina- 
tor of the vector meson propagator an imaginary part to 
deal with this fact. 

The result of Eq. l|11.7fl is valid above the two-Goldstonc- 
boson threshold. At P 2 — 4M?, the 'regular part', Eq. I|11.8fl . 
vanishes, and remains zero below the threshold, since as 
we remarked above the integral I then has no regular part 
and is completely 'infrared singular'. The two represen- 
tations for the infrared singular part, valid for different 
ranges of the parameter P 2 , may be 'joined together' at 
the threshold singularity. A similar thing happened in the 
last chapter for the two representations of the infrared 
singular part of the scalar loop integral Iv</>- 



11.2 Application to the self-energy 

The major problem in finding the 'soft part' of the ampli- 
tude of Fig. |31 has been solved in the last paragraph. For 
the full expression, we need to add some vertex structure 
from the local effective Lagrangian. We choose to work 
with the interaction Lagrangian of Eq. (|2.11|l and refrain 
from constructing interaction terms with a higher num- 
ber of derivatives, though not all momenta in the present 
problem can be considered as 'soft'. Since the coupling- 
constant Gy may be measured from p-meson decay, where 
the Goldstone bosons are also not of soft momentum, this 
can be seen as a valid approximation. Applying the usual 
Feynman rules, we obtain 



2 F 4 

d d k [k^P v -k v P^){PPk a 



^ j^jj^v,pa x jabc jbad 



P a k p ) 



(27T)* (k*-Ml){{k-PY-Ml) 



(11.9) 



Before further evaluating this, we have to discuss the power 
counting. The vertices are both of the order 0(p), since 
only one momentum in the product k ■ P is a small mo- 
mentum in the sense of the power counting scheme. Re- 
membering the discussion of the last paragraph, we want 
the amplitude to be of 'chiral order' e?+l + l — 2 = d. We 
will see that the infrared regularized amplitude will indeed 
respect this power counting. Using the tensor integral of 
App. El we get 



. — 
V zed r>fiv,pa 

2F 4 ] 



(11.10) 



l_(^ + i(4M|-P 2 )7 



Here I is the scalar loop integral of Eq. (jl 1 . 1|1 . and we 
defined 

p^u.pa = g^ppvpa _ g pa p u p p _ ffl ^ v y (11.11) 

The infrared regularized amplitude is obtained from 111.9fl 
by simply letting / — > Iir. The infrared part Jir, was of 
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0(p d ~ 2 ). In order to check that the terms of 0(p d ~ 2 ) can- 
cel in the soft part of H11.9JI . it is easiest to use that the 
first term in the chiral expansion of Iir, is also the first 
term of the series for / so ft, which was given in Sect. HTTP 

2 

^soft — ~j^I<i> + • • • 

Inserting this in {—i)^v'm, that is the infrared part of 
Eq. Ijll .9(1 . it is clearly seen that the infrared part of the 
amplitude is indeed of order 0(p d ), as required by low- 
energy power counting. 



11.3 Contributions to the Vector Meson Mass 

First we introduce some notation. We define 

1 = l^' pa = ^(g^g 1 " 7 - g pa g vp ). 
Furthermore, we write 

p = p^v,po 

see Eq <|ll.llfl ■ It is easy to calculate 

11 = 1, 1 P = P , 
P 1 = P , P P = 2P 2 P , 

where the multiplication works as e.g. 

1 a ■ P Q / 3 'P <T — V>P a 

The tensor field propagator may then be written 



(11.12) 



D 



M 2 



21 



M 2 - pa J' 



(11.13) 



while its inverse (in the sense of the above multiplication) 
is 

D-'-i(^l-ip). 

The one-particle irreducible self-energy amplitude may be 
parametrized as 



E = ^41-*P 

2 4 



(11.15) 



where A and B are scalar functions of P 2 and the meson 
masses. 

The procedure is now standard: Summing over the 
number of self-energy insertions, we find that the full prop- 
agator 

D fuU = D + D(-i)SD + . . . 

is given by 



D 



full 



(D^+tE)- 1 
( 



M 2 (l-A) 



21 



V 



M 2 (^)-P 2 



(11.16) 
\ 

/ 



We have to look for the poles of this expression. Since A 
is a small perturbation of 0(p 2 ), the only pole will be at 



M 2 



1 



B 



(11.17) 



with My iP h the physical mass of the vector meson. Before 
we use this formula to compute the contribution of Fig. [3] 
to the vector meson mass, let us make a very rough esti- 
mate of the expected size of the contribution. The most 
general effective Lagrangian for the tensor field contains 
a term 

c(W^W^ X +), 

yielding, among other terms, a contact term contribution 
of 0(p 2 ), which gives rise to a shift of the propagator pole: 



Mi 



Mi 



8cMi 



Since the coupling constant c is not known, for the pur- 
pose of our estimate we make a naturalness assumption 
concerning this coupling, and set c = 1, which gives us 
a value of 100 MeV for the mass shift. If power counting 
is a consistent perturbative scheme here, we would ex- 
pect for an 0(p 4 ) correction a number of size of roughly 
{Ml/Ml) (100 MeV) ~ 3 MeV (for the pion contribu- 
tion). Now let us compare this estimate with the (infrared 
regularized) amplitude corresponding to Fig.yj It will con- 
tribute to B, defined above, with 



R 6G v 

LJV — T~ 

V pi 

1 



4Mi 



M 2 



12 



'IR 



1447T 2 



(d l M^-(l+Ad 1 )M^)\, (11.18) 



giving a mass shift of 1.2 MeV, which is really only a 
small correction, and also of the size expected by the (very 
rough) estimate made above. If we had used the full (real 
part of) the integral I, we would get a result that is com- 
parable to a correction of 0{p 2 ) (of course, there are such 
terms of 0(p 2 ) in the full integral, i.e. the power-counting 
violating terms). We conclude that the main effect of the 
graph of Fig. U (at the physical pion mass) is due to the 
imaginary part of this diagram, associated with the width 
of the vector meson propagator. 



12 Chiral extrapolation of the rho meson 
mass 

In this section, we analyse the quark mass dependence of 
the p-meson mass and related topics. This is not entirely 
new, see e.g. Refs. [31,32], but we do not want to rely 
on any model or the assumption of 'dominating' contribu- 
tions to the p self-energy. In fact, there are many different 
contributions to the self-energy of the vector mesons, and 
only a few of the corresponding LECs are known from phe- 
nomenology. Of course, one could resort to models like the 
massive Yang-Mills approach or the extended NJL model 
to estimate these parameters (as it is done e.g. in the 
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work of Bijnens and collaborators [25]), but our goal is 
more modest. We resort to parameterizing the pion mass 
dependence of M p and fix the combinations of LECs from 
existing lattice data [33]. This allows e.g. to analyze the 
value of M p in the chiral limit. 

First, let us discuss the many different contributions 
to the vector meson mass. We restrict ourselves to terms 
at most quadratic in the quark masses. The first type of 
contribution stems from tree diagrams with quark mass 
insertions, i.e. operators ~ x+ or ~ X+j like e.g. 



(W-W X+ > ,(W-W)( X - 



(W-W 



(12.1) 



The LECs accompanying such explicit symmetry breaking 
terms are in general difficult to determine, as it is well 
known from the analysis of the nucleon mass in chiral 
perturbation theory, see e.g. [20,34,35]. Such tree graphs 
lead to the following vector meson mass terms: 1 



M$ ee 



h Ml + k 2 Ml 



(12.2) 



with fci (fe) a combination of dimension two (four) LECs. 
There is also a tree graph without quark mass insertion, 
it corresponds to the vector meson mass in the chiral lim- 
ited, denoted as My in what follows. Next, we consider the 
various one-loop graphs. Tadpole diagrams with an inser- 
tions of the second order effective chiral Lagrangian have 
also to be considered, some of the pertinent structures are 

(W • W x+) , (W • W u a u a ) , {W^W^g^ u a u p ) 

Note that in addition to the symmetry breakers of the 
type given in Eq. I|12.1(l . kinetic terms ~ d^d^cj) from 
the second order effective Lagrangian also contribute, thus 
increasing the number of LECs to be determined. In the 
comparable case of the nucleon mass, these can be deter- 
mined to good accuracy form the analysis of pion-nucleon 
scattering in the low energy regime. The total contribu- 
tion of the tadpoles to the vector meson mass takes the 
form 

M Wpole = fe3M 4 ln f^ ( (123) 

with fc 3 another combination of dimension two LECs. The 
sunrise diagram (cf. Fig. ^l) starts to contribute at order 
p 3 because there are one derivative vertices of the form 

{e^ pa Wp V V a W ap u a ) , 

A famous example of such a vertex is the uipw coupling, 
which is generated in meson field theory from the Wess- 
Zumino-Witten term, see e.g. [11,12]. It was e.g. consid- 
ered in the analysis of [32] as one of what these authors 
call 'dominating contributions'. Since there are various of 
such VV4> couplings, we write the sunrise contribution to 
the vector meson mass as 



My DrK 



I M 2 



My 



(12.4) 



To avoid notational clutter, we absorb all prefactors like 
1/F 2 etc. in the coefficients fci. 



which is again reminiscent of the leading non-analytic con- 
tribution to the nucleon mass. The ellipsis denotes analytic 
terms ~ M% and higher order contributions. Finally, we 
have to consider the self-energy graph considered in the 
preceding section. It leads only to a fourth order contri- 
bution of the form 



(12.5) 



To be specific, we consider now the pion mass expan- 
sion of the p-meson mass, i.e. we set My — M p and 
M$ = Mi; m the above formulae. Including only the non- 
analytic terms from the fourth order, it takes the form 



M p = M° 



c x Ml 



c 2 Ml- 



■c 3 M*ln{j^)+0(Mt), 
(12.6) 

where M° is the mass in the chiral limit, and the Ci (i = 
1,2,3) are combinations of coupling constants as discussed 
before. In the absence of a detailed phcnomenological anal- 
ysis of these couplings, we will use the CP-PACS data [33] 
for the p-meson mass as a function of the pion (average 
light quark) mass to determine the parameters M^, ci, c 2 
and C3. We only employ lattice data with M 2 < 0.5 GeV 2 . 
In fit 1, we fit these parameters by demanding that the 
physical p-mass is obtained for M^ = 140 MeV. For fits 2 
and 3, however, this restriction is lifted. In these fits, we 
input the chiral limit mass. Throughout, the fits are sub- 
jected to the further restriction that one obtains natural 
values for the combinations of LECs, that is we enforce 
\ci\ < 3. The corresponding fit parameters (obtained by 
least-square fits) are collected in Tab. ^ 
The corresponding curves are shown in Fig. To get a 
better handle on the theoretical uncertainty, we also allow 
the fits to stay within the theoretical uncertainty of the 
lowest point at M 2 = 0.1 GeV 2 , as shown in Fig. [5J If we 
insist again on naturalness of the coupling constants, we 
can bound the p-mass in the chiral limit by 



650 MeV < M° p < 800 MeV 



(12.7) 



These results are similar to what was found in the pioneer- 
ing work in Ref. [32], but they are less model-dependent. 
The range for is also consistent with the numbers de- 
rived by Bijnens and collaborators in their study of vector 
mesons in chiral perturbation theory [25]. It would be in- 
teresting to extend these studies in two directions, first to 
include also more recent lattice data and second to try to 
give more stringent limits on the combinations of LECs 
by incorporating more phenomenological constraints. 



Table 1. Fit parameters. * denotes an input quantity. 





Fit 1 Fit 2 Fit 3 


M p u [GeV] 
ci [GeV" 1 ] 
c 2 [GeV" 2 ] 
c 3 [GeV" 3 ] 


0.776 0.650* 0.800* 
-0.662 2.200 -1.215 

1.291 -1.934 1.915 
-1.723 1.572 -2.367 
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0.2, 0.3 

M„ 2 [Ge V s ] 

Fig. 4. The rho meson mass as a function of the light quark 
mass, M% ~ (m u + m<j). The solid (dot-dashed) line(s) refers 
to fit 1 (2,3) as described in the text. The lattice data are from 
CP-PACS [33]. The diamond denotes the physical rho mass. 




0.2, 0.3 

K [GeV 4 ] 

Fig. 5. The rho meson mass as a function of the pion mass: 
Theoretical uncertainty as described in the text. For further 
notations, see Fig. 2] 



The quark mass expansion of the p-mass EQ. (|12.f)|) al- 
lows one to deduce the corresponidng cr-term, 



A dM p M2 dM p 



* dMl ' 



(12. 



with rh the average light quark mass. From the numbers 
collected in Tabled we nn d 



1.9M2<<7^<1.5Mj 



(12.9) 



This shows again that the rho as a massive particle has a 
very different quark mass expansion than the pion, where 
ov ~ M% [1]. In magnitude, the rho cr-term is similar to 
the pion-nucleon one, a^N — 45 MeV. 



perturbation theory. The pertinent results of this investi- 
gation can be summarized as follows: 

1) The most economic way to deal with vector mesons 
in chiral perturbation theory is to utilize the anti- 
symmetric tensor field formulation as stressed in [13]. 
When vector mesons appear in tree graphs only, cal- 
culations are straightforward as summarized in Sect. [3 
and App.fA] Of course, other formulations like the vec- 
tor field approach can also be used, see App. IBICI 

2) When vector mesons appear in loops, the appearance 
of the large mass scale complicates the power count- 
ing, as discussed in Sect. |31 and Sect. 01 in essence, 
loop diagrams pick up large contributions when the 
loop momentum is close to the vector meson mass. To 
the contrary, the contribution from the soft poles (mo- 
menta of the order of the pion mass) that leads to the 
interesting chiral terms of the low-energy EFT (chiral 
logs and alike) obeys power counting. We have briefly 
summarized the method proposed in [19] to extract the 
'soft pole' contribution from one-loop integrals. 

3) The standard case of infrared regularization [20] , where 
the heavy particle line is conserved in the (one-loop) 
graphs is recapitulated in Sect. For these cases a 
very elegant splitting of a Feynman parameter integral 
allows to unambigouosly separate the infrared singular 
from the regular part, cf. Ea. i|5.5[l . 

4) In the case of spin-1 fields, new classes of self-energy 
graphs appear. The case for lines with small exter- 
nal momenta but a vector meson line appearing inside 
the diagram in analyzed in Sect. and the singularity 
structure of the corresponding integrals is discussed in 
Sect. [7| In Sect. [8] the infrared singular part for such 
types of integrals is explicitly constructed, cf. Eq. (|8.4() . 
As explicit examples, the Goldstone boson self-energy 
and the triangle diagram are worked out in Sect. ED and 
Sect. El respectively. 

5) A different type of one-loop graphs appears in the vec- 
tor meson self-energy, where only light particles (Gold- 
stone bosons) run in the loop. This is discussed in 
detail in Sect. 1111 where the corresponding infrared 
singular part is extracted, see Eq. i|11.7|) . and the con- 
tribution to the vector meson mass is worked out. We 
briefly discuss the problems related to the imaginary 
part of such type of diagrams. 

6) As an application, we consider the pion mass depen- 
dence of the p- meson mass in Sect. El We show that 
there are many contributions with unknown LECs, still 
one is able to derive a compact formula for M p (M 7T ), 
see Eq. <|12.6[) . We analyze existing lattice data [33] 
and conclude that the p-meson mass in the chiral limit 
is bounded between 650 and 800 MeV. We have also 
discussed the irp sigma term. 



13 Summary and outlook 

In this paper, we have considered chiral perturbation the- 
ory in the presence of vector and axial- vector mesons (spin- 
1) fields and presented an extension of the infrared reg- 
ularization scheme originally developed for baryon chiral 



The methods outlined here can be applied to many in- 
teresting problems, for example one could systematically 
analyze vector meson effects on Goldstone properties like 
form factors or polarizabilities or extend these consider- 
ations to systems including baryons (for a first step see 
e.g. [23]). ' 
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A Tensor field approach 

In this appendix, we briefly discuss the representation of 
spin-1 particles in terms of antisymmetric tensor fields, 
following closely Appendix A of [13]. 

As an antisymmetric tensor field, W pu has six degrees 
of freedom, whereas a massive vector field only has three. 
Loosely speaking, there are two spin-1 fields 'hidden' in a 
general antisymmetric tensor field approach (correspond- 
ing to a reducible representation of the rotation group). 
To make this clear, we decompose the tensor field W pv (in 
momentum space) 

W pv = (W-PWP) pv + (PWP) pv = W^ + W^, (A.l) 
where the matrix P is the projector 



(A.2) 



and p p is the momentum four-vector associated with the 
tensor field. Because of the projector property of P, we 
have 

P"W£ = 0, 

giving 4 conditions for six degrees of freedom, but one 
condition is redundant due to the antisymmetry property 
of W^,,. So we are left with 6 — 3 = 3 degrees of freedom 
for W T , and therefore also for W L . 

Inserting the above decomposition in the general form 
of an action principle for antisymmetric tensor fields [13], 



S w = J d 4 x{(a - 2b)8 p W pu d p W pv + bd p W pv d p W» v 
+cW pu W^} , (A.3) 
with arbitrary parameters a, b, and c ^ 0, and using 

wj^wp = 0, 

which is easily verified by a direct calculation, we see that 
the action splits in two terms: 



->w 



where 



S Wl = J d 4 x{ad p W^d p W p » + cW£,Wr), 
S Wt = J d 4 x(bd"W^d p W^ + cW p T v W pv ). 
Therefore, the path integral can also be factored: 

J[dW]e i§w = J[dW L ]e iSw t J ' [dW T ]e iSw T . (A.5) 



In the rest frame, this decomposition corresponds to 

3 3 

[dw\ = [dw L ][dw T ] = n^oiin^-]- 

i=l i<j 

Following [13], we choose b — 0, so that the second path 
integral becomes an unimportant constant (from the view- 
point of the classical action, Wt becomes a non-propagating 
field). For a = 0, in contrast, Wl would be 'frozen' in this 
way. Choosing, furthermore, 



1 Ml 

a = — t:> c = — 1~- 

2 4 



(A.6) 



and dropping the letter L, the massive vector field is de- 
scribed by 

J [dW]e lSw = J [dW}exp^i J d A x {-^d ll W pv d p W pv 



M 2 

+ ^W pv W^)}- (A.7) 

Without invalidating the above argument, we can also add 
a coupling linear in W, of the form 



with some external antisymmetric current J^ v . Since we 
do not use couplings quadratic in the tensor field in this 
work, we will not attempt to extend the given argument to 
Lagrangians including such more complicated interaction 
terms. 



B Vector field approach 

One can now ask how much the results obtained from cal- 
culations similar to the ones in the main body of the text 
depend on the description in terms of an antisymmetric 
tensor field. Of course a dependence of that kind is not 
wanted for physical observables! 

What happens if one uses a more conventional vector 
field approach for the description of the vector mesons? 
The Lagrangian for a massive vector field V p is well known: 

c k v m = -hy^vn + l -Ml(V p V») (B.l) 



(A. 4) where is the field strength tensor associated with V p , 



v pv = v p v v -v v v^ 

The corresponding vector field propagator (in momentum 
space) reads 



G pu (k) = H) 



k 2 - M 2 + it ' 



(B.2) 



If one now tries to write down interaction terms in analogy 
to the ones given before, it is seen that they are all 0(p 3 ), 
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giving a resonance exchange diagram of 0(p 6 ) in contrast 
to the 0(p ) result derived with the antisymmetric tensor 
field description. In particular, the interaction terms are 



pint 
l-V 



fv 

2V2 



igv 
2V2 



{K,u v ]Vn- (B.3) 



Here fy and gy are new coupling constants, and the minus 
sign is a pure convention. This interaction looks very much 
like Eq. (|2.11(l . but note that the tensor V^„ contains, from 
its definition, an additional derivative giving the interac- 
tion the order 0(p 3 ) instead of 0(p 2 ). The form factor 
contribution derived with this interaction is of higher or- 
der than the result of the last section and does not agree 
with experiment. If we had started with a vector field de- 
scription for vector mesons, we might have concluded that 
the vector meson contribution is less important than sug- 
gested, for example, by a dispersive analysis. This problem 
was already solved in [15], we simply repeat here some of 
salient ingredients in a slightly different way. 

The mathematical relation between the two variants 
of the theory was worked out in [15] by imposing e.g. the 
large momentum transfer constraints on the pion vector 
form factor. In App. [U] we argue that the Lagrangians 



and 



nkin 



4M, 



16 My 



t jni. 



(B.4) 



(B.5) 



indeed give equivalent theories on the level of the path 
integrals, when J M „ are some traceless hermitian sources. 
This equivalence is achieved by simply integrating out the 
field W M „, which leads to the appearance of contact terms 
quadratic in the source. The source terms we use are 

V = ~V2(F V F+ + »Gv[u M , u u \) 

(compare Eq. (|2.11|) 1 and, assuming the validity of the 
path-integral argument for the equivalence, we draw the 
conclusion that equivalence leads to 



Fy 
My' 



gv 



Gy_ 

My ' 



(B.6) 



in perfect agreement with Ref. [15]. Furthermore, since 
the V uv - couplings are of 0(p 3 ), we conclude that vector 
meson exchange to lowest order 0{p A ) can be represented 
by the contact terms 



1 



16M V v " 



G 2 



! " AM? ^' 



SM V 



Using the definition of the objects u u and F+ v and stan- 
dard trace relations, the contact terms can be written as 

+L V P 9 + LY0P10 + HYxPn , 



with Pi the structures of the fourth order meson Lagrangian 
[2] and furthermore 



LX = —QL\ 



lX = 


Gy 
8My ' 


^2 — 


2LY, 




FyGy 




Fy 


2My ' 


AMy 



TTV _ Fy 
H ll 



8My 



C Duality transformation 

The action for an antisymmetric tensor field W UL , is 



S w = I atx C w (C.l) 



Note that now stands for the field Wjj v , having three 
degrees of freedom, as discussed in App. ^] Now let us 
modify this action by a term containing a vector field V^: 

1 M 2 (A 

(C.2) 

We have also added a source term linear in W uv . The equa- 
tions of motion (e.o.m.) following from this Lagrangian are 



+ " W = 2Mj/^ (C - 3) 



(C.4) 



Substituting V n from the second equation into the first 
one, the latter becomes the 'old' e.o.m. following from 
Lw — \ J ■ W alone. So on the classical level, the addi- 
tional squared term does not change anything. 

Considering now path integrals, we deduce by a Gaus- 
sian integration over V that 



[dW] / [dV]e 



if d A xLv,w 

(C.5) 

where the '~' means 'up to a constant factor'. It must 
be noted here that the vector field V should be treated 
like W, that is using a constraint to eliminate one degree 
of freedom. On the classical level, the second e.o.m. says 
that V has the same number of degrees of freedom as W, 
i.e. three d.o.f. Contributions of the path integral over V 
deviating from this e.o.m are exponentially damped like 
e~ x because of the form of the squared term in Ly^w- 
Using a constraint on V, we will get just another constant 
prefactor as long as this constraint is linear in the field. 

If we now recklessly interchange the order of integra- 
tions on the right-hand-side of Eq. (|C.5(I . define 

Vu V = d^V u - d v V^, 

and do the lU-integration, we get 



[dW] 
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J[dV][dW]e X p{i Jd 4 x (^M V W^W^ - i J^W^ 
J [dU]expjz Jd 4 x (-hy lu ,v>» + ^-V^ 



+ 4M v JtwV 16M 2 ^ 

(C.6) 

This corresponds to a conventional vector field Lagrangian, 
together with contact terms of the form 



1 



t nil/ 



16M 2 ^ 

If the tensor field W is given in the matrix notation of 
Eq. I]2.6H , which implies that one has to take the (flavor) 
trace of every term in the Lagrangian L\y, the source term 
linear in the field W^ v can be taken as 



4V2 



If J^y is traceless and hermitian (which is the case for 
the couplings we consider in Ea. (|2.11[l L we produce the 
contact terms 



1 



1 



■■ r v ) V2 (jFf) V2 (rf) ~ lour"" 



which justifies the form of the contact terms given in the 
preceding appendix. 

A path-integral approach to show the duality of the 
antisymmetric tensor field and the vector field description 
has also been presented in [36], but the approach used 
here is shorter. In principle, we have not done much more 
than a Legendre transformation to change variables from 
dW to V. The PF-field has been integrated out, leaving 
as a 'souvenir' only the contact terms. As in the foregoing 
Appendix, we do not attempt to include also more com- 
plicated coupling terms in such a transformation, with the 
excuse that in this work the contact terms are only needed 
for the linear interaction of Eq. I|2.11|l . 



D Loop integrals 

In this Appendix we treat the reduction of general loop 
integrals to linear combinations of scalar loop integrals. 
The scalar loop integrals and ly have been defined in 
Eas. (|9.5[) and (|9.6() . respectively, while Iv$ was defined in 
Eq. (GlU). Writing q 2 = (q 2 



M 2 ) 



My , one derives 



d d c 



(2tt)<* (q 2 -M 2 )((q-p) 2 -M 2 ) 



1$ + MyI V<t> (p). 



(D.l) 



Moreover, by Lorentz invariance, we must have 
d d q q» 



(2n) d (q 2 M 2 ){{q - p) 2 - M 2 ) ^ (D ' 2) 

because there is no other four-vector than p^ available 
here. The scalar S can be found by contracting with p^: 

a = J_ f A q 2 + p 2 - (g - pf 

2p 2 J (2ttY {q 2 ~ M 2 )((q~ p) 2 - M 2 ) 

= ^(l4>-Iv + (M v -M 2 +p 2 )I V4> ). (D.3) 

Using both results together, we can compute 
d d q (p ■ q)q» 



{2-kY (q 2 -M v )((q-p) 2 ~M 2 ) 



= -^~2 (a><t>I<t> + aylv + av<plv<p), 

where the coefficients are given by 

a^ = M 2 -M 2 + 3p 2 , 

a v = M 2 - M v - p 2 , a V4> = , 

These results so far are already sufficient to arrive at the 
decomposition of Eq. (|9.3|) for Is . Another useful result is 
also derived using Lorentz invariance: 

d d q 



(D.4) 



(D.5) 



(27T) d (q 2 -M 2 )((q-p) 2 -M 2 ) 



gTA{p) 



p»p v 

9 

V 



B(p) 



(D.6) 



The coefficients can be found by contracting with g^ v and 
Pn, using the above results and g l _ LV g tJ " v = d : 

(d - l)A(p) = ^ + i(4M| - p 2 )I^(p) (D.7) 
(d- l)B(p) = (~ - l)l 4 + (±p 2 - Mljl^(p). 

Here the integral I^{p) is obtained from Iy^ip) by sub- 
stituting M^ for My- 

The last result we have to derive is the integral I^{p, fc), 
which we only need for on-shcll kinematics , p 2 = k 2 = 
M|. Using the algebraic decomposition 

4(p ■ k)q 2 - 4(p ■ q)(k ■ q) 
(q 2 -My)(q 2 -2q-p)(q 2 ~2q.k) 

4(j3 • k)My - My 



(q 2 -M 2 )(q 2 -2q-p)(q 2 ~2q-k) 



4(p -k)-q 2 - M, 



(q 2 -2q-p)(q 2 -2q-k) 
1 1 



M 2 



q 2 — 2q ■ p 



2q ■ k 



My 



+ 



(q 2 -M 2 )(q 2 ~2q-p) 
M 2 

( q 2 -M 2 )(q 2 ~2q-k) ' 
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the results obtained in this appendix enable us to calculate 
I T A (p,k) = 

d d q {p + k- 2q) T ((p ■ k)q 2 - (p ■ q)(k ■ q)) 
{2ttY (q 2 - M 2 ){{q - p) 2 - M 2 )((q - k) 2 - M 2 ) 



= (p + k) T (dv^Iv^iPi k) + d v< pl v< p 
+d^ >< pl ( p < p(k - p) + d^I^ + dyly) 



(D. 



The coefficients are given by 

{Ml -M v + (p- fc))(4(p • k)M v - My) 



Ml + (p-k) 



Ad 



v</> 



Ad d 



4cL 



2M^M V - M v _ M v - 4(p ■ k)My 
Ml Ml + (p-k) ' 

M| + (p • fe) 4(p • fc)Af £ - My 



d- 1 

1 Ml 



Ml + (p-k) 



My 



Xr, 4d v ='—\- 1. 

d-i Mr y mi 



(D.9) 



We remind the reader that this is valid only for p 2 = k 2 
M?. For this case, we have Iy<p(k) = Iv<j>[p) = Iv<f> ■ 
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